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1 Abstract 

We construct a theory which has the same particle content as softly broken minimal supersym- 
metric QED (MSQED) and is free of quadratic divergences up to two loops. Also this theory is 
completely gauge invariant. It appears that MSQED is not at all the only theory without these 
quadratic divergences. This proves that there exist non supersymmetric theories in which there 
are no quadratic divergences up to two loops. 
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2 Introduction 



In this paper we discuss an extension of QED where all quadratic divergences are absent up to 
two loop order. Of course the softly broken supersymmetric extension of QED (MSQED) hase 
this nice property. The Lagrangian of this theory can be constructed with the help of superfields 
(see e.g. ref. [T]). In this way one can find the Feynman rules for MSQED. These Feynman rules 
are given in appendix 1. 

For supersymmetric theories one can prove ([2], 0) that all Green's functions are free of 
quadratic divergences (QD). This is considered to be one of the great merits of supersymmetric 
theories, i.e. supersymmetric theories solve the naturalness problem. It is then natural to ask 
whether the specific supersymmetric theory one has constructed is the only theory without QD. 
It is this question that we address in the case of MSQED. 

In [5], [S] and [7J similar questions are addressed. Veltman [3] showed that in the Standard 
Model the QD vanish at 1-loop order if one imposes one constraint on the coupling constants. 

Inami et al. [5] investigate a model with one Majorana fermion and a number of real spin-0 
fields. They demand that no QD occur up to 2-loop order. They find that for a theory with 2 
or less spin-0 fields this condition uniquely determines supersymmetry, i.e. the supersymmetric 
theory is the only theory without QD. 

Deshpande et al. [B] investigate several models and impose that all QD have to vanish at 
1-loop order and that the equations one gets for the coupling constants should be invariant under 
the renormalization group equation. They find, in all the cases they consider, that with these 
constraints one always gets the supersymmetric theory. Because they demand that their 1-loop 
equations are invariant under the renormalization group their calculations are actually of higher 
loop order, i.e. they demand that QD vanish at all orders in perturbation theory. 

Jack et al. [7J investigate theories where, if the QD vanish at 1-loop order, the QD at 2-loop 
automatically vanish. They give an example of such a theory which is not supersymmetric. 

In [5] and [5] the absence of QD up to some loop order means the theory under consideration 
is necessarily supersymmetric. In this paper we show that this is not the case for all theories. Of 
course we do not claim that there exist non supersymmetric theories without QD at all orders 
in perturbation theory, i.e. we do not claim that we can solve the naturalness problem without 
introducing supersymmetry. It is just interesting by itself how many orders one needs before one 
gets to a supersymmetric theory (if one gets there at all). As is shown in this article, in which we 
limit ourselves to 2-loop calculations, two loops is not enough to get to the supersymmetric theory. 
Also for this reason we don't demand that our equations are invariant under the renormalization 
group, (as they do in e.g. [S] and [13].) If we would do so our calculation would effectively be of 
more than 2-loop order. To keep a clear sight on what happens at every order in perturbation 
theory we explicitly calculate amplitudes at every order. Of course practical reasons have disabled 
us to carry on with our analysis at 3-loop order and further 

This paper is organized as follows. We start by postulating the particle content of our theory, 
which is the same as MSQED. Then we start constructing the vertices one by one. We will 
only introduce new vertices when it is absolutely necessary to keep everything free of QD. When 
introducing a new vertex we will keep the form general. Then we will fix it as much as possible 
by insisting that there are no QD. This process will be done up to 2-loop order. 

After this construction up to 2-loop order we will see how much of the theory is actually fixed. 
It will appear that, to keep all Green's functions free of QD up to 2-loop order, the form of the 
vertices are fixed, however there is still a lot of freedom to choose the value of coupling constants. 
At this point we will also check that the theory we have built is gauge invariant. 

3 Building the Theory 

We wish to build a theory with the same particle content as MSQED. This means we have one 
spin-i fermion with charge e and mass m e (electron e), one spin-1 boson with mass (photon 
7), two scalar particles with charge e and mass mi and niR (left and right selectron &l and ejj) 
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and one spin-^ Majorana fermion with charge and mass m,y (photino 7). Because the photino 
is a Majorana fermion special Fcynman rules have to be used. A convenient formulation of these 
Feynman rules is given by Denner et al. in [8]. These rules are explained in appendix 1. Note 
that we work in the Feynman/Lorentz gauge. In section [7] we will come back to the point of 
gauge invariance of our theory. 

Now we want to build a minimal theory with these particles which is gauge invariant, unitary 
and free of QD. Let's start by introducing the electron photon interaction we know from ordinary 
QED: 




(1) 



If we only have this vertex we have QED with three extra particles which do not interact. 
We know then that the theory is gauge invariant and unitary, but not free of QD. In the vacuum 
polarization at 1-loop order we encounter a QD. Only one diagram, , contributes to this 
vacuum polarization, and its QD part is easily calculated: 



n 
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QD 



(2tt)4 
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■ IE 



ie 



(2) 



In the last step we kept only the QD part. How we compute the QD part of a diagram is explained 
in appendix 2. We see that only the longitudinal part of the photon propagator contains a QD, 
the transversal part is free of QD because it is protected by the gauge invariance of the theory. 

Here we see the need of another vertex with a photon. The only 4 possibilities are depicted 
below. 







Now for the moment, to keep our theory as simple as possible let's only introduce the first two. 
The general form of these vertices is: 



E^ = Ap^ + Bq 1 " 



(3) 



where p M is the momentum of the incoming selectron, q^ is the momentum of the outgoing 
selectron and A and B are constants. All momenta are counted from left to right. Of course A 
and B could depend on p 2 and q 2 , however to keep our theory as simple as possible we will not 
introduce form factors in the vertices. The exact form of the two vertices can be fixed with the 
Ward Takahashi identity (see e.g. [3]) for the 3- vertex at tree level: 



-E^q-p) 



p* — m* L ' ' 11 q 2 — m" L p" 
Substituting the general expression ([3]) in (j4]) we find: 

A = -ie, B = - 

so that both vertices become: 



ie 



(4) 



(5) 




-ie{p + qY 




Mp + qT (6) 
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Now, to see if we still have unitarity consider the process ~^ 77- The Ward identity 

(also see e.g. [§]) states that if we take one of the photons longitudinally polarized and off shell, 
but the other particles on shell, the matrix element should be zero. There are two diagrams: 




All indicated momenta are again flowing from left to right. Taking photon 1 longitudinally 
polarized [e\ — qi) and off shell a simple calculation gives: 

Mi + M 2 = -2ie 2 qi ■ £ 2 • (8) 

This means we need another diagram to satisfy the Ward identity. Now, it can easily be seen 
that none of the possible 3-vertices one can construct is going to help here. The simplest thing 
we can do then is introduce a new 4- vertex. It is clear that the only 4- vertex that will help is: 




This new vertex gives one extra diagram and it is easy to see from ([HJ that it has to be of the 
form A g 1 *" , with A = 2ie 2 . Of course the same reasoning can be given for the process ene-R — * 77 
and we get two new 4- vertices: 




Now we go back to the photon propagator. Because of the new vertices we get 4 new diagrams, 
of which the QD parts are readily calculated with the method outlined in appendix 2: 



Note again that only the longitudinal parts of all diagrams get QD. We see now that all QD 
cancel in the sum Py V + . . . + P^ v . So with these 4 new vertices we solved the problems of the 
longitudinal part of the photon propagator. 

Now consider the left selectron propagator. With the vertices we have introduced up to now 
we have two contributing diagrams: 



L 




L L 



The QD do not cancel, so we need a new vertex again. It is easy to see that the only 3-vertices 
which will help here are vertices where a left selectron, an electron and a photino meet. There 
are four of these: 




Because the Lagrangian should be Hermitian there are two relations between these four vertices: 




(12) 



Now because of these new vertices we get one new diagram: 




Let's call the first of our new vertices B>l(p, q): 



B L (p,q) = 




where p is the momentum in the electron line, flowing from left to right and q is the momentum 
in the photino line, also flowing from left to right. Now the general expression for this vertex is 
a Fierz decomposition: 

B L (p, q) = SI + V al a + T a(i a^ + A a7 V + Pi 5 (15) 
where a a P is defined as 

<^ = 5(7V-VV). (i6) 
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Since we do not want any form factors in our fundamental vertices, to keep the theory as simple 
as possible, S and P should be constants and for V a , A a and T a p we can write down the general 
expressions 

V a = Vi p a + v 2 q a 
A a = ai p a + a 2 q a 

T af 3 = t 1 paqf}+t 2 p l 3qa+h£af35 1 p 5 q 1 (17) 

where also v\, v 2 , fli, a 2 , ti, t 2 and t% are constants. Notice we haven't included a term propor- 
tional to 

Qa/3 in T-cx.fi since 

g af 3 a a0 = . (18) 

Now P3 is given by 

U 1 f jijrv ( U I 1 1\ i{1>-y+ m e) 

P 3 = ~7^4 / d lTr \ B L(p-l,-l)-? pro TT^' 

(27r) 4 J I (p-iy-mj + te 

<- 1 )/^(P-'.-07»|tt^) ( 19 ) 

Note that we don't have to worry about the special Feynman rules for Majorana fermions here, 
since there is only one electron line the fermion flow can be fixed along the electron arrow 
unambiguously. Now the general expression (fT5|) should be substituted in (fT9|) . This involves a 
lot of algebra and can best be done by computer. It appears we find even quartic divergences 
(i.e. proportional to A 4 ). The coefficients in front of these terms have to vanish, from which we 
find two equations: 

-3|ii-i 2 | 2 -12|i 3 | 2 = 0^h=t 2 , < 3 = 
4\ Vl ~ v 2 \ 2 + 4\ ai - a 2 \ 2 = 0=>v 1= v 2 , a 1 =a 2 (20) 

When we substitute these results in the coefficients for the QD parts proportional to p 2 and (p 2 ) 2 
we find one more equation (the coefficient for (p 2 ) 2 A 2 is automatically zero when we substitute 

my- 

M 2 + k| 2 = o (21) 

which means that v\ = v 2 — a\ — a 2 = 0. Note that because t\ = t 2 the whole third term in (|15p 
vanishes because of the anti symmetry of . 

So our general expression l|15p is now reduced to 

B L = SI + P7 5 . (22) 

For S and P we get one constraint because the quadratically divergent part with no factor (p 2 ) n 
in front should cancel the QD in P\ and P 2 : 

-3e 2 +4(|Sf + \P\ 2 ) = . (23) 

Of course one can also consider the right selectron propagator and impose the absence of QD. 
Then a relation similar to (|2"3"|) is found. So finally our 4 new vertices become: 



G 



with the constraints 



<-» S L I + P Ll 5 



-S L I + PlI 5 



SrI + P fl7 5 



-S R I + Prh 5 



(24) 



3e 2 + 4(|5 £ | 2 + |Pi| 2 ) = 



-3e 2 + 4(|S fl | 2 + |P R | 2 ) 







(25) 



To get one more constraint from one loop processes we can consider — ► Sr. There is only 
one diagram for this process: 



Pi= M^}^ ^ (S* l Sr + PIPr)^A\ (26) 

so that we also get 

S* l Sr + P* l Pr = 0. (27) 

With the vertices we have up to now one can also calculate the QD in the electron and 
photino propagator. It appears that these have no QD, because of their fermionic character. 
Also diagrams with 3 or more external legs, which we haven't considered up to now, will have no 
QD, which can be verified by simple power counting. 

So now we have constructed a theory free of QD up to 1-loop order. We still have freedom to 
choose our coupling constants Sl, Pl, Sr and Pr (8 real parameters against 4 real equations) 
and we haven't even seen the need for introducing the selectron 4- vertices that occur in MSQED. 
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4 QD at 2-Loop 



Now let's consider our theory at 2-loop order and impose the absence of QD. First consider the 
left selectron propagator. There are 29 1PI diagrams: 






Note that for these 2-loops diagrams the fact that the photino is a Majorana fermion is important. 
Otherwise diagrams 14, 15, 24 and 25 wouldn't have been possible at all. It is also for these 
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diagrams that the special Feynman rules of the Majorana fermion become important because the 
arrows of the electron lines clash. 

Because of the large number of diagrams we have automatized the computation of the QD in 
the propagators. The program LOOPS. frm, written in FORM [TO] can do the calculation of the QD 
for us. This code generates all diagrams at two loop order (also the 1PR ones) and keeps only 
the QD parts. The QD parts are written in terms of the standard integrals A and B: 



(2tt) 4 J * l\ + is (h - l 2 ) 2 +ie I 2 + is 



ie) 2 l\ + is 



(29) 



Notice that we omitted the masses, since these are unimportant for the QD part of the integral, 
integrals A (or B) with different masses in the propagators can just be added as far as their QD 
part is concerned. Also the external momentum p can be put to zero in the denominators, since 
they arc also unimportant for the QD parts. 

Running the program in the case of the selectron propagator gives us the coefficients in front 
of A and B. It appears that the coefficient of B already vanishes if we substitute the constraints 
we found at one loop order. The A-part gives a new constraint, its coefficient is: 

= llie 4 -2ie 2 \S L \ 2 ~2ie 2 \P L \ 2 -l0i\S L \ 4 -m\P L \ 4 -l6i\S L \ 2 \P L \ 2 + 



HSi) 2 p 2 + 2iS 2 (p* L ) 2 - m\s L \ 2 \s R \ 2 - m\p L \ 2 \p R \ 2 



-2i\S L \ 2 \P R \ 2 - 2i\P L \ 2 \S R \ 2 + 2iS* L P L S R P R + 2iS L P* L S R P R + 
-6iS* L P L S R P R - &iS L P* L S R P R (30) 

We could now proceed with the other propagators at 2-loop order, however, first there is a 
problem. Consider the three constraints we found at one loop order again: 

3e 2 -4|S L | 2 -4|P L | 2 = 
Ze 2 -A\S R \ 2 -±\P R \ 2 = 

S L S* R + P L P R = 0. (31) 



From the second equation we see that 



I^| 2 = 7e 2 -I^| 2 - (32) 
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From the last one it follows that 

Pl = ■ (33) 

Substituting both relations in the first equation in (|3ip gives: 

\Sl\ = \Pr\ • (34) 

From (l3"3l it then follows that 



|5 fl | = |Pi|. (35) 
If we use all these results in the 2-loop condition (|30|) this condition reduces to: 

^e 4 = . (36) 



Apparently our constraints are incompatible. This forces us to introduce new vertices. The only 
possible 3- vertices we can still introduce are: 



With the Ward Takahashi identities these vertices can be fixed like we fixed ©, they will have 
exactly the same expression as ©. But this means we will ruin the photon propagator at one 
loop order again! 

So we have to introduce new 4-vertices. Since the problem occurs with the selectrons it seems 
plausible to introduce selectron 4-vertices. These are also the simplest 4-vertices we can introduce 
because of the scalar character of the selectrons. So we will introduce: 




-iVr 




-iV R 



(37) 




-iV L R 



(38) 



Because of these new vertices the left and right selectron propagator at 1-loop order also 
change. It is easy to see that our new 1-loop constraints become: 



3e 2 
3e 2 



MSl\ 2 



4|fY 



V L + V LR = 
4\S R \ 2 - 4\P R \ 2 + V r + V lr = 
S L S* R + P L P R = 



(39) 



Nothing changes in the electron, photon and photino propagator at one loop order. 

The new constraint for the selectron propagator at two loop order can also calculated with 
LOOPS. frm again. We get: 







llie 4 -2ie 2 \S L \ 2 -2ie 2 \P L \ 



2i{SiYPi + 2iSi(Pi) 
-2i\S L \ 2 \PR\ 2 -2i\P L 



10i|S L | 4 -HM|P L | 4 -16i|S L | 2 |P L | 2 

m\p L \ 2 \p R \ 2 + 

2iS* L P L S R P R + 2iS L P* L S R P R + 



m\s L \ 2 \s R \ 2 



-6iS* L P L S R P R - d>iS L P* L S R P R - 2ie 2 V L - 2ie 2 V LR 



2t\S L \ 2 V L - 
+ iV, 



2i\Pl?V l + 2i\S R \ 2 V LR + 2i\P R \ 2 V LR 



LR 



The QD part of the right selectron propagator and the process 
also be calculated. From these we get: 



(40) 

e R at two loop order can 



= 



2ie'\P R \ 2 - m\S R 



10z|Pfl| 4 

wi\s R \ 2 \s L \ 2 ~ m\p R \ 2 \p L \ 2 



m\s R \ 2 \p R \ 



Hie 4 - 2ie 2 \S R \ 2 - 2ie 2 
2i(S R ) 2 P 2 +2iS 2 R (P^ 2 
-2i\S R \ 2 \P L \ 2 - 2i\P R \ 2 \S L \ 2 + 2iS* R P R S* L P L + 2iS R P* R S L P* L 
-6iS* R P R S L Pl - &iS R P* R S* L P L - 2ie 2 V R - 2ie 2 V LR + 
2i\S R \ 2 V R + 2i\P R \ 2 V R 



2i\S L \ 2 V LR 



2i\P L \ 2 V LR 



^V 2 



iV? 



LR 



(41) 
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= -2ie 2 SiS R -2ie 2 P* L P R -%i\S L \ 2 P* L P R -m\S L \ 2 S* L S R + 

2iS l {p* l ) 2 s r - %i\p l \ 2 s* l s r - %i\p r \ 2 s* l s r ~ m\s R \ 2 s* L s R + 

2iS* L S* R P 2 R + 2i(S* L ) 2 P L P R - m\P L \ 2 P* L P R - %i\S R \ 2 P* L P R + 

2iP* L S 2 R P R - m\P R \ 2 P* L P R + 2iS* L S R V LR + 2iP* L P R V LR (42) 

Of course (|4Tj) is just (j40f with L and R swapped. 

Also the QD parts of the electron, photon and photino propagator can be calculated with 
LOOPS. frm. These appear to be zero already, independent of our choice of coupling constants. 



5 1 Loop Conclusions 

Using our constraints from 1 loop ([SU]) we can eliminate 4 of our 11 free (real) parameters. Let's 
define: 



S L = 


\Sl\ 


e *t>L 


Pl = 


\Pl\ 




s R = 


\Sr\ 


e i4>R 


Pr = 


\Pr\ 





(43) 

Now we use ((3^|) to express everything in the 3 phases <j> R , ipi and tp R , the magnitude \Pr\ and 
the 4- vertex constants Vl, Vr and Vlr- 



+ IpL — i>R + It + 27T7Z 



Here a is defined as 



and the inequalities 



\Sl\ 2 = a\P R \ 2 

\Pl\ 2 = a(^e 2 -\PR\ 2 + \v R + \v L R 
\Sr\ 2 = \e 2 -\PR\ 2 + \v R + \v L R (44) 



i 



\V R + \V L R 



V l + V L r > -3e 2 

Vr + V L r > -3e 2 (46) 

should be satisfied in order that there exists a solution at all. 

We see there is still a lot of freedom left if one only wants a theory free of QD up to 1-loop 
order. 



6 2 Loop Conclusions 

Now we can use the 1 loop relations (|44l) in the 2 loop constraints (|4"U|) . (|41[) and (|4"2"|) . The easiest 
way is to use the 1-loop relations in (|4"2"|) first. Then this equation simplifies to 

(SlSr P L P R ) (l i) ((P£) 2 - (^2) 2 ) = . (47) 

We can satisfy this equation in 3 ways: 

1. S L S R = P L P R 

2. a = 1 

3. P* L = ±S* L 
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6.1 Case 1 



In this case (|40| and (|41jl simplify to two quadratic equations in Vl, Vr and Vlr, it appears that 
\Pr\ drops out. We can use these two equations to express Vl and Vr in terms of Vlr- This 
gives 4 solutions. We are left with 4 real degrees of freedom: 



i>L, tf>R, \Pr\, V L r (48) 



The 4 solutions are: 



ip L = i'L + -7T + 7r(m + n) 

lPr = ip R - -7T + 7r(m - ra) 

\Sl\ 2 = |Pfl.| 2 

1^1* = W - \? - \Pr\ 2 + ± J v HOe 4 + 44e2Vz.fi - 2V? R 



and 



Vl = Vr = lle 2 + V L Ji± \/ 11 0e 4 + 44e 2 V Lfl -2V i 2 ii (49) 



1 

— 7r + 7r(m + n) 
1 

— 7r + 7r(m — n) 





= IpL ~ 


VR 


= i>R - 


S L \ 2 


= a\Pi 


Pl\ 2 


- °G 


Sr} 2 


= Ie> 
3 


V l 


19 , 

= y e 


Vr 


19 , 



3 e 2 - |Pfi| 2 + -V L * T 397e4 + 226e2y LR ~ 11V L 



§ r 1 

2 + \V LR T ^V 397e4 + 226e2 ^« ~ 11V tfl 



'r\ 2 + ~V L r T f^\/ 397e4 + 226e2v &* - n ^fl 



with 



_ 28e 2 + 8V Lfl ± ^397e 4 + 226e 2 V L fl - 11V^ H 
" ~ 28e 2 + 8V Lfl T V397e 4 + 226e 2 V Lfl - HV L 2 fl 

6.2 Case 2 



(50) 
(51) 



In this case ()40[) and (|4ip simplify to the same equation, where also |Pr| drops out. This equation 
has two solutions in terms of Vlr- We are left with 5 real degrees of freedom: 



VR, ipL, ipR, \Pr\, V L r (52) 



The two solutions are: 



Vl = VR + i'L - ipR + 7T + 2-7rn 
\S L \ 2 = \Pr? 



\Pl\ 2 = \Sr\ 2 = 7 -e 2 - \Pr\ 2 + \v LR ± \ ^HOe 4 + 44e 2 V Lfl - 2V L 2 fl 



Vi = V R = lle 2 + V Li? ±^/ll0e 4 + 44 e 2 V L fl-2V i 2 fl (53) 

We see that the two solutions in this case are like the first two solutions of case 1 (jM]), just 
slightly more general (Here ipR is also free.). 
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6.3 Case 3 



In this case we also get two quadratic equations in Vl, Vr and Vlr, which also have four solutions. 
We are left with 3 real degrees of freedom: 



and the four solutions are: 



ipL, tpR, V LR (54) 



ip L = ipL + 

ifR — tjjR — 7r + 7rm 



LR 



and 



V L = V R = lie/ + V L r ± J n0ei + MeW L R-2VZ R (55) 



f>L = i>L + Km 

<PR — l/lR ~ K + 

I^I 2 = I^lI 2 = + \v LR ± y a v 183e4 + 132e2 ^ - &v Ir 

\Sr\ 2 = \Pr\ 2 = 7 -e 2 + \v LR T y a V 183e4 + 132e2 ^ - 6V Zr 

V l = 4e 2 + V LR ±jsj 183e 4 + 132eW LR - Wl R 



Vr = 4e 2 + V L r T \ y/ 183e 4 + 132e 2 V LR - 6Vg R (56) 

The first two solutions are plotted in figure [TJ In this figure the crosses indicate the values 
for the MSQED coupling constants, they lie on the curves of the (-) solution (|55|) . 

The final conclusion is that, in the case of MSQED, MSQED is not at all the only theory 
free of QD up to two loops. With the demand that up to two loops all QD vanish one can fix 
the form of the vertices in a minimal extension of QED. For the values of the coupling constants 
there is still a lot of freedom left. 



7 Gauge Invar iance 

While building our theory we have fixed some of our vertices by using the Ward or Ward Takahashi 
identities. We chose our vertices such that the processes under consideration satisfied these 
identities. This does not mean however that our complete theory is gauge invariant. To see that 
our complete theory is gauge invariant for general Sl, Pl, Sr, Pr, Vl, Vr and Vlr we will 
demonstrate this at the level of the Lagrangian density. From our Feynman rules the Lagrangian 
density can be read off. We find: 

C = C K + Li , (57) 
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|SJ 2 /e 2 (+) 
V L /e 2 (+) 

|SJ 2 /e 2 (-) 
V L /e 2 (-) 



Figure 1: The solutions (55). \S L \ 2 /e 2 = \P L \ 2 /e 2 = \S R \ 2 /e 2 = \P R \ 2 /e 2 And V L /e 2 = V R /e 2 
are plotted as functions of Vlr/c 2 . Both the solution with (+) and (-) are shown. The MSQED 
solution is indicated by the crosses. 



with 

C K = -^F^F^-^{-i^d^ + m^^-il{a + m 2 L )i L -i R (a+m 2 R )i R + 
-Ipe (-Z7 M 9 M + m e ) tp e , 

Ci = -eA^jf/tpe + 

-ieA* [e^df^L - (d^A e L ) - ieA 1 [e R d^e.R - ( d^e R \ e R ^j + 
e 2 A 2 e\e L + e 2 A 2 e R e R + 

-i e L ip e (-S* L + P£ 7 5 ) V 7 " i A^i ( s l + PlI 5 ) ipe + 
-i e R 4> e (-S* R + P^ 7 5 ) - i 4^7 (Sb + p fl7 5 ) V'e + 

e[e L e[e L - -Vr e R e R e R e R - V LR e f L e L e R e R , (58) 

where of course 

F» v ee d^A v - d v A^ . (59) 

Now it is easy to check that indeed this Lagrangian density is invariant under the gauge 
transformation: 

A fl (x) -» A"(x)+d^A(x) 
ip e {x) — » exp{-iek(x))ip e [x) 
&l{x) — > exp (— ieA(x)) &l 

sr{x) — > exp (— ieA(x)) §r (60) 

(Strictly speaking it is not this Lagrangian density that is invariant, but the action S — J d A x C, 
the Lagrangian density itself is invariant up to terms which are total derivatives.) 
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A Appendix 1: The Feynman rules for MSQED 



In the three figures below the Feynman rules for MSQED are given in the Feynman/Lorentz 
gauge. The straight line with arrow indicates the electron, the combination of a straight line and 
a wiggly line indicates the photino, the straight lines without arrow with a label L or R indicate 
the left or right selectron and the wiggly line indicates the photon. Here the photino is a Majorana 
fcrmion, which makes it impossible to fix a fermion flow in each diagram unambiguously. To be 
able to work with Feynman rules we follow Denner et al. [5]. They state that for any fermion loop 
one should choose an arbitrary orientation for the fermion flow. One should count the momenta 
along this chosen fermion flow and the vertices T change as: 

T — ► CT T C _1 . (61) 

Here T denotes transposition and C is the charge-conjugation matrix. The expression on the 
right hand side is equal to 

CT^C" 1 = rjT (62) 



where rj is given by: 



. i for r = i,n 5 ,yV (63) 



The fermion flow in all diagrams will be indicated by an arrow next to the fermion line. 

Figure 2: The Feynman rules for the bare propagators in MSQED. The indicated momenta flow 
from left to right. The arrows next to the fermion lines represent the chosen fermion flow. 

p z — m*-\-te 



~ 2 'i i - 



p 2 — m~ +ie 



p 2 — m 2 R -\-iE 



) 2 +ie 
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Figure 3: The Feynman rules for the 3- vertices in MSQED. The indicated momenta flow from 
left to right. The arrows next to the fermion lines represent the chosen fermion flow. Notice that 
the electron-selectron-photino vertices where the fermion flow is opposite to the electron arrow 
are not explicitly given, these vertices have the same expression as their counterparts where the 
fermion flow is along the electron arrow. 





-ie(p + qY 




-ie(p + qY 



- 75(1 + 7 5 ) 




- -7f(I-7 5 ) 



Mi-7 5 ) 



V2 




- -7f( I + 7 5 ) 



1G 



Figure 4: The Feynman rules for the 4-vertices in MSQED. 
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B Appendix 2: Computing QD parts 



To see how a loop integral diverges one has to introduce a regulator, in our case a cut-off mo- 
mentum. This cut-off can be a strict cut-off at momentum A, which means that no momentum 
occurring in a bare propagator may exceed A. In our case it is however more convenient to work 
with a smoother cut-off momentum ([11]. [12] ) introduced by 



1 



1% + m? J a- 2 



da e - a V 2 E+ m2 ) (64) 



where Ie is a Euclidean momentum. This means that in all our loop integrals we have to do a 
Wick rotation first, to get a Euclidean form. 

The loop integrals (1 and 2 loop) for propagator diagrams which we have to know will in 
general depend on the momentum p flowing through the propagator and several masses. However, 
if one is only interested in the worst divergent part of a loop integral these can in principal all be 
set to zero. We shall do this for the momentum p, but not for the masses, to avoid introducing 
infrared singularities. On can however take all masses equal. 

Now consider the 1-loop integral 

In = ! dH m -j ; — . (65) 



ie) n 

After the Wick rotation one finds: 

In = i(~l) n I d% - 2 - „. (66) 



(1% + m 2 ) 1 

Using the regularization scheme introduced above this becomes: 

I n = nr 2 (-l) n / dU!... da n j e -(a 1 + ...+a„)m 2 ^ 

J a- 2 (a% + . . . + a n ) 

From this expression it is easy to see the worst divergent parts: 

h = -ITT 2 A 2 

7 2 = 2in 2 In A 

In>2 = ^ 2 (-l) n - -1 TT^ZI, (68) 

(n — 2)(n — 1) m zn 4 

indeed verifying what one expects from simple power counting. 
The only 2-loop integral which we will encounter is: 

/"111 

1=1 d i lid A l2 75 5 — n TT^5 5 ~ 75 5 ~ • (69) 

J If — m 1 + ie [li — 12) — m z + ie t| — m z + ie 

After Wick rotating, using the regularization and performing the momentum integrals this be- 
comes 

I = ir 4 d ai da 2 da 3 e -(ai+a 2 +a 3 )m 2 _ ^ 7Q ) 



(ai«2 + aic*3 + a2as) 2 

The worst divergent part (in this case quadratically divergent, as we will see below) can be 
obtained by putting the exponential to 1. Then our 2-loop integral becomes 



(|) A*. 171) 



indeed showing a QD. 
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